ON NONEXISTENCE OF BARAS-GOLDSTEIN TYPE 
FOR HIGHER-ORDER PARABOLIC EQUATIONS 
WITH SINGULAR POTENTIALS 



V.A. GALAKTIONOV AND IV. KAMOTSKI 

Abstract. The celebrated result by Baras and Goldstein (1984) established that the 
heat equation with singular inverse square potential in a smooth bounded domain C 
R N , N > 3, such that Oefi, 

u t = Au + j^u in ftx(0,T), u\ gQ =0, 
in the supercritical range 

C > CHardy(l) = 

does not have a solution for any nontrivial L 1 initial data uq(x) > in O, or for a positive 
measure. Namely it was proved that a regular approximation of a possible solution by a 
sequence {u n (x,t)} of classical solutions of uniformly parabolic equations with bounded 
truncated potentials given by 

V(x) = j^^ V n (x) = min{^, n) (n > 1) 

diverges, and, asn-^ oo, 

u n (x,t) -> +oo in f2x(0,T). 

In the present paper, we reveal the connection of this "very singular" evolution with a 
spectrum of some "limiting" operator. The proposed approach allows us to consider more 
general higher-order operators (for which Hardy's inequalities were known since Rellich, 
1954) and initial data that are not necessarily positive. In particular it is established 
that, under some natural hypothesis, the divergence result is valid for any 2mth-order 
parabolic equation with singular potential 

u t — ~(—A) m u + | a ,| C 2m u in f2x(0,T), where c > ch(«i), m > 1, 

with zero Dirichlet conditions on d£l and for a wide class of initial data. In particular, 
typically, the divergence holds for any data satisfying 

uq(x) is continuous at x — and uq(0) > 0. 

Similar nonexistence (i.e., divergence as e — * 0) results are also derived for time- 
dependent potentials e~~ 2m q(j, j^r) and nonlinear reaction terms ^J"^^ with p > 1. 
Applications to other, linear and semilinear, Schrodinger and wave PDEs are discussed. 
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1. Introduction: Baras-Goldstein result and extensions 



The paper treats the questions of regular approximations of higher-order parabolic 
equations and the corresponding elliptic operators with singular unbounded potentials. 
There are two main generating and key ingredients of such a theory: 

(i) The origin of such research can be attributed to Hardy's inequalities (1919-20) for 
such symmetric operators in the L 2 -space, which established the parameter supercritical 
range, where the operators are not semibounded and hence do not admit Friedrichs' 
classic self-adjoint extensions (1935). It is also crucial that, at the same time, in view of 
coinciding deficiency indices, there exist infinitely many other self-adjoint extensions with 
discrete spectra and L 2 -eigenfunctions. 

(ii) It was later the discovery of Baras-Goldstein (1984), that, in the corresponding 
non-stationary parabolic flows, regular approximations of nonnegative solutions uniformly 
diverge, i.e., the corresponding semigroup gives infinite values for any such nontrivial data. 

It turned out later on that, in general singular linear and nonlinear PDE theory, these 
are important problems concerning existence of extended semigroups of proper (blow-up) 
solutions that are obtained by regular approximations. For a number of parabolic, hy- 
perbolic, and other evolution equations of mathematical physics, such semigroups can be 
essentially discontinuous in any suitable or admissible metric. Actually, this means that 
a proper solution blows up at some t = T > 0, ceases to exist as a bounded solution 
for t > T, so a special framework of extended semigroup theory via regularization should 
be put in charge instead of the classic one. In particular, proper setting of basic prob- 
lems such as the standard Cauchy or initial-boundary value ones represents a difficult 
subject. Concerning related elliptic operators with singular potentials, it is then key to 
understand, which of their spectral properties admit a proper approximation by a family 
of "regularized" operators with truncated singularities. It turns out that, even for classic 
Laplacian-type operators with inverse square potentials in the supercritical Hardy range, 
standard self-adjoint extensions have nothing to do with the actual operator that occurs 
in the regularization limit. It is revealed that this actual "limiting" operator has positive 
eigenvalues and we demonstrate that they are responsible for the effective "nonexistence" 
of the limit of a solution. This point of view allows us not to rely on the Maximum 
Principle (which was a standard tool in previous work related to Baras-Goldstein type 
problems) and to consider much more general operators for which the Maximum Principle 
is not available. In particular, we demonstrate it considering higher-order elliptic opera- 
tors. Such 2mth-order operators are common in PDE theory and mathematical physics, 
and the first derivation of Hardy's inequality for m = 2 goes back to Rellich (1954). In 
what follows, we explain these aspects in greater detail treating also other PDEs. 

Thus, one of the first key results of modern theory of singular elliptic operators and 
related extended discontinuous semigroups of blow-up solutions was obtained by Baras 
and Goldstein [I]. 

1.1. Baras-Goldstein (1984): nonexistence for singular potential in linear heat 
equations. Let Q C M. N be a smooth bounded domain, where N > 3 and G fl. In 
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Baras-Goldstein the authors considered the heat equation with the inverse square 
potential and the zero Dirichlet boundary condition, 

(1.1) u t = Au + j^u in Qx(0,T), u(x, t) = on dQ x (0, T), 
and nonnegative initial data 

(1.2) u{x, 0) = u {x) > in Q, 

where uq G L 1 (fi), or, is a positive measure. 

Their key nonexistence result is as follows: in the supercritical Hardy range 

(1.3) c>c H = (^) 2 

the problem (11. ip . (11. 2p does not have a solution^. Namely, this means that the regu- 
lar approximation of a possible solution by the sequence {u n (x, t)} of bounded classical 
solutions satisfying (11.11) with the truncated bounded potentials obtained on replacement 

(1.4) V(x) = ™ h-> V n (x) = min {™, n}, with n > 0, 
diverges. More precisely, as n — > oo, 

(1.5) Un(x,t) -> +oo in fix(0,T). 

On the contrary, for c < ch in (II. ip . this sequence has a finite limit that corresponds to 
existence of a unique solution of the Cauchy problem (11.11) . (II. 2p . 

The results and ideas of the pioneering paper [TJ generated a new direction of nonex- 
istence/existence theory for linear and nonlinear PDEs. We refer to [TTl [T8| I2T] for the 
study of linear parabolic equations and to [161 122] devoted to quasilinear reaction-diffusion 
PDEs. These questions are reflected in the books [25] and Ch. 11]. 

1.2. Main extensions to higher-order PDEs without positivity assumptions on 
data. As our main extended model, we consider the 2mth-order parabolic equation with 
a singular potential, 

(1.6) it* = B ?/ = -(-A) m u + d&jr it in Q x E + , where m > 1, N > 2m. 
For definiteness, we take zero Dirichlet conditions on the smooth boundary 

(1.7) « =§* = ... = f^¥ = on fflxl +) 
where i/ is the unit outward normal to dQ, and initial data 

(1.8) u(x,0) = u (x) E L 2 (Q). 

Actually, as one can expect, the nonexistence result is purely associated with the strong 
singularity at x = and does not essentially depend on boundary conditions (though 
their self-adjoint nature is a convenient assumption for calculus applied). 

As in Baras-Goldstein [1] , our crucial assumption is that the potential in (I1.6P belongs 
to the supercritical range, i.e., the constant c is large enough: 

(1.9) c> c H (m). 



The question on nonexistence was posed to the authors by H. Brezis and J.-L. Lions, [TJ p. 122]. 
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Here cn(m) for dimensions N > 2m is classic Hardy's best constant of multiplicative 
inequalities involving the potential in (jl.fip . This goes back to Hardy (1919) for m — 1, 
and Rellich (1954) for m = 2; see [T3j and |29j for further references and full history. The 
Hardy constant is given by 



;i.io) chM 



B 2 B i ...B m for m even, 

B 3 B 5 ...B m ch(1) for m odd. 



Here = j^ ( jV ~ 2fc )(^ r+2fc ~ 4 ) j 2 f or £ = l ; 2,...,m, and ch(1) is as in (jl.3P ; see a simple 
derivation and a list of references in [8] , and also [5j [26] with a number of applications for 
m = 1. 

Similar to [1] and as usual in extended semigroup theory [HIE], we construct a proper 
solution of the problem (ll.6l) - (ll.9p using regular approximations. For convenience, in- 
stead of (11.41) (clearly, this does not affect the final results), we perform an analytic 
approximation of the potential by replacing 

(1.11) V{x) = ^jk i-> V e (x) = £2m . + %| 2m , with e > 0. 

By {u £ (x,t)} we denote the sequence of classical bounded (for t > 0) solutions of this 
regularized initial- value problem for the parabolic PDE with the potentials (II. lip . 

(1.12) u £ : ut = B £ u = -(-A) m u + £2m ^ x[2m u in SlxK + (m > 1, N > 2m). 

We then take the same Dirichlet boundary conditions 111- 7ft and same initial data (11.81) . 

Passing to the limit as e —>■ is performed in Section [3], which is based on a more general 
approach in Section [2] to 2mth-order parabolic equations with arbitrary potentials. In 
particular, we show that the phenomenon: 

(1.13) {u £ (x,t)} diverges in L 2 (f2) as e — >• 0, 

is a generic (robust) property of such approximations of singular parabolic problems. 
As a simpler counterpart, it follows that for c < ch the limit of {u £ } always exists, 

(1.14) u £ (x, t) — > u(x, t) as e^O, 

for any data uq G L 2 (Q), where u(-,t) G L 2 (Q) for all t > 0. This means existence of a 
(unique) solution. 

The nonexistence result (11.131) actually means that, in the supercritical range c > ch, 

(1.15) \u(x, t)\ = oo 

can be treated as a "proper solution" of the original problem (ll.6l) - (ll.8p for any L 2 
data satisfying some extra conditions, e.g., the regularity and positivity at the origin 
(or, more generally, non-orthogonality to a positive lineal, see below). Then, as we have 
mentioned at the beginning of Introduction, (I1.15P means that the extended semigroup of 
such solutions obtained by regular approximations is discontinuous at t — for sufficiently 
arbitrary initial data uq. Note that the behaviour of the "approximating" family {u £ (x, t)} 
as e — ► to get (11.151) can be extremely oscillatory For instance, for data uq(x) that are 



oscillatory at x = 0, we present in Section H] an example of radial solutions with non- 
uniform oscillatory blow-up, where u E {x,t) has both limits ±00 along some subsequences, 
and the same holds even at the origin x = 0. 

In the present paper, we study the case where the singularity of the potential is concen- 
trated at the single internal point x = G Q, and this indeed affects our final nonexistence 
conclusions. We do not consider the case of singularities on the boundary, which can lead 
to other regularization spectral properties and hence different and more difficult nonex- 
istence criteria. In the elliptic case, there exist related extended Hardy's inequalities for 
smooth domains Q corresponding to potentials with inverse square singularity on the 
boundary dQ, which were introduced in [I]. See also [TO], where further extensions via 
perturbations of such potentials were presented, and also [TJJ [T5] for delicate nonlinear 
counterparts. 

1.3. Layout of the paper. Sections EH6] are occupied with various aspects of nonexis- 
tence analysis of parabolic models such as (11.61) and their time-dependent and nonlinear 
extensions. To show further application, in Section we discuss some rather simple 
corollaries of our analysis being applied to Schrodinger and hyperbolic PDEs such as 

i u t — — (— A) m w + i^prr u and u tt = -(-A) m u + u, 

together with their semilinear counterparts. 

2. General nonexistence theorem 

2.1. Divergence theorem for a general potential. In this section, we treat the nonex- 
istence in a more general setting. We consider the Cauchy-Dirichlet problem for the 
poly-harmonic equation with a general regularized potential, 

(2.1) u £ : u t = B £ u = -(-A) m u + ^q(^)u in flxl + (m > 1). 

Here the potential q e (x) = ?(f ) depends on the parameter e > 0, for instance, in 
a manner similar to that in (j!.12p . where q(y) is now an arbitrary smooth function in 
~R N decaying at infinity. Therefore, (12 .ip plays a role of a regular approximation of the 
parabolic equation with the singular potential V(x) (V(0) = 00 and |V(a;)| < 00 for 
i^O) such that, uniformly on compact subsets in \ {0} and sufficiently fast, 

(2.2) q £ (x) = ^ 9 (f ) -> V(x) as e - 0+ (e.g., V(x) = ^ as in <^M) . 
We assume for u E (x,t) zero Dirichlet conditions on the smooth boundary, 

(2.3) « e = ff = - = f^ = on dQxR + , 
where v is the unit outward normal to dQ, and prescribe initial data 

(2.4) u s (x,0) =u 0e (x) e L 2 (Q) 
that, in general, also depend on e. 
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In what follows, C and c £ denote various positive constants which exact values are of 
no importance. Here, C is independent of e, while constants c e have at most rational 
dependence on e, i.e., 

p — p 

e < c e < e , 

for some positive P which can be arbitrarily large. In this calculus, we easily write 
c £ c £ > c £ . Such constants are not of importance and are negligible while dealing with 
exponential factors such as e 1/,£ , e 1//£ , or e 1 ^ m , to be treated more carefully. 

First of all, using the above calculus, since the equation is linear, we may always assume 

(2.5) c e < ||woe|U 2 < c £ for e > small, 

i.e., uoe does not get large as e — * 0, which is a natural and not restrictive assumption. 

Our main hypothesis on the potential is as follows: 
Hypothesis (P): The limiting operator. 

(2.6) B = —(—A y ) m + q(y)I in R N 

has M > 1 positive eigenvalues with exponentially decaying eigenfunctions, 

and, for some its eigenfunction U n (y), 

(2.7) \(u 0s (z),U n (z)}\ > c £ for all small e > 0. 
In the following particular cases, < \2.7\j is replaced by: 

(i) if uo £ (ex) — > 5q > as e — > (e.g., uq is independent of e, Uq(x) is continuous at 
x = and Uq(0) = 5q > 0), there exists U n (y) that has nonzero mean in M. N ; and 

(ii) if uoe{x) is supported in some ball B 5s (0), with ^ — ■> and j u 0e (x) > c e; t/iere 
exists an eigenfunction U n (y) of <\2.6§ that does not vanish at the origin y = 0. 

Our main nonexistence result is a corollary of the following estimate: 

Theorem 2.1. Let m>l, ( 1 1 . S\\ . and Hypothesis (P) hold. Then the sequence {u E (x,t)} 
of classical solutions of the approximating problem <\2. ip - pL|D diverges in L 2 (Q) as e — > : 

(2.8) ||w £ (x, £)||l2(q) > e^ 2 ™' — > cxo for any fixed t > 0. 

2.2. Proof of Theorem 12.11 Fix a sufficiently small £ > 0. Consider the corresponding 
eigenvalue problem 

(2.9) B £ ip = Xi/j, H 2m {Q) n # m (^)- 

By cr(B e ) = {A^, j > 0} and {ipj, j > 0} we denote the corresponding spectrum and 
orthonormal, complete, and closed in L 2 (Q) eigenfunctions subset, [2]. Clearly, 

(2.10) u e {x, t) = J2u> 0) c| ^( x ) eAlt ' with c ) = ( M o e , ^ £ ), 

(2.11) so that |k(-,*)||i 2(n) = E (i > ) \ctfe 2X j l . 
We will need the following statement: 



Lemma 2.1. Let Aj, Uj for j = 1,2, ...M, be positive eigenvalues and corresponding 
normalized eigenf unctions of (12.61) . Then, for any small enough e > 0, intervals 

l £ 2m " e ) £ 2m ' " £ J 

contain eigenvalues of the problem (12. 9p . Moreover, the function Uj = Uj{^) is an 
"approximate" eigenfunction of (Iff. ,9ft in the following sense: 

(2-12) ll^ £ -E(i)«^ll^(Q) <e"T, 

where E(i) stands for the summation with respect to j such that X e j belongs to the above 
interval, and E(i) a ] = c e> a j being constants (depending on e). 

Let us return to the proof of the theorem. Due to the assumptions of the theorem, we 
fix the necessary eigenfunction U n of the operator ( 12. 6p that provides us with the estimate 

(2.13) \(u e,U n )\ >c e . 

On the other hand, due to the above lemma, we see that there is at least one eigenfunction 
ipf of (12.181) (see below) with the corresponding eigenvalue Af in the interval 

Xfe [^-Ce-T^ + Ce-T]. 

In view of (12.121) we have, 

(2.14) \(uo £ ,U n )\-e-T\\u 0£ \\ < | £(i) #>| < (£<i)K«fc.Vf>r)*, 
and consequently, (12.51) and (I2.13P immediately implies that 

E(l) \ c l\ 2 > c s- 

Then obviously we have 

\\u £ (x,t)\\l m = E >o) I^W > E(i) Itf |W > cj^ E(i) I4I 2 > c £ e^. 
This concludes the proof. □ 

2.3. On generic asymptotic behaviour as e — > 0. The following sharp "pointwise" 
estimate of the solution sequence {u £ } holds: if (12. 7p is valid for j = 0, then, for arbitrarily 
small fixed t > 0, 

(2.15) e- A «V(s,f)-Co#)^0 as e^O in L 2 (fi) (ifj £ (x) ~ e~% U (f) ). 

2.4. A weaker blow-up hypothesis. The blow-up result ( 12. 8 p remains valid under the 
weaker condition on initial data: (I2.5P holds and 

(2.16) \{u 0£ (x),Uj(-))\ > e~~^ for some positive constant c* < Cj, 

where the constant Cj gives the rate of exponential decay of the eigenfunction U n , i.e., 

(2.17) WM^Wh^) < C i with c i >0 - 

Assumption (I2.16P is obviously satisfied in the pioneering original paper [Tj for j = 0, see 
Proposition 13.11 (iii). 
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2.5. Proof of Lemma 12.11 Consider eigenvalues and eigenfunctions of operator (12.61) 
Aj, Uj, j = 1,2, M. These eigenfunctions decay exponentially at infinity, see (I2.17p . 
Then functions satisfy the equation 

(2.18) B £ ^(f) =^A^(f) in a 

These do not satisfy boundary conditions on dQ. To fix this, consider the functions 

V j {x) = U j (£) X (x), 

where % is smooth function which is equal to 1 in a neighborhood of the origin and is 
equal to zero in some neighborhood of the boundary dQ. Then 

(2.19) B e Vj = AjVj + Vj in fi, Vj G H 2m (n) n # m (fi), and 

(2.20) \\V 3 \\ Hm{n) <e-T. 

Then application of the "Lemma on approximate eigenfunction" (see [27]) and observation 
that the difference between Uj and Vj can be estimated via the right-hand side of (12.201) 
deliver the result of the lemma. □ 

Lemma [2. II can be strengthened (though, in this form, we are not going to use it): 

Proposition 2.1. Let Aj and Uj are eigenvalues and eigenfunctions of the problem ( 12.61) . 
Then eigenvalues and eigenfunctions of the eigenvalue problem satisfy: for any fixed 
l<j<M, ase^O, 



where the summation in ) t a ^ es place over k such that Ak = Aj, and £( K .)(a|) 2 



C F . 



3. Problem for ( 11.121) : scaling, limit problem, and limit operator 

We now explain in greater detail application of the above general results to our main 
model equation (11.61) with the approximation (11.121) . 

3.1. Scaling. Thus, we perform in (II. 12ft the natural scaling of the independent variables 

(3.1) x = ey, t = e 2m T. 

This reduces (11.121) to the equation with the regular potential, 

(3.2) Ut = ^ iU = -^^ u + _^_ u in Q eXR+i 

which is posed in the domain 

(3.3) n £ = {yeR N : ey G fi} H aQe = 0). 
The initial data are now 

(3.4) u 0e (x) i ^ u 0e (ey) in Q £ . 
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3.2. Limit problem. The above parabolic problems fl3.2l) - fl3.4p are posed in a family of 
expanding domains {fi e } that embrace the whole R N as e — > 0. Therefore, it is natural 
to consider the limit problem, which is the Cauchy problem for 

(3.5) Vr = Bl v = -(-A) m v + TT ^v in R N xR + . 

3.3. Spectral properties of limit operator. The limit parabolic problem for (13.51) 
introduces the limit operator 

(3.6) Bl = -(-A) m + TT ^J in L 2 (R N ), where c> c H (m). 
Obviously, f)3.6p is semibounded, 

(3.7) B 1 <cl in C?{R N ), 

so it admits Friedrichs' self-adjoint extension (denoted again by Bi) [2J with the domain 
H 2m (R N ) in view of the embedding [231 P- 54] 

(3-8) j < C J RN \D m w\ 2 in H m (R N ). 

For convenience, we present a detailed description of the necessary spectral properties: 

Proposition 3.1. (i) The spectrum of the operator (\3.6§ with domain H 2m (R N ) comprises 
the continuous and the discrete ones, 

(3.9) <r(Bi) = <r c (Bi) U <r p (Bi) = {A < 0} U {A = A,- > 0, j = 0, 1,2, ...}, with 

(3.10) c> A > A 1 > A 2 > ... > A k > A fc+1 > ... 0, 

where each eigenvalue repeated as many times as its finite multiplicity. 
(ii) Each eigenfunction Uj{y) is exponentially decaying at infinity, 

,l/2m, I 

(3.11) \Uj(y)\ < Aje- amA i M as y -»• oo, 

where Aj and a m are positive constants. 

(hi) The first eigenfunction Uo(y) is radially symmetric. For m = 1, it is positive, 

(3.12) U (y) > in R N (m = 1). 
For any m > 2, Uo(y) has infinitely many sign changes. 

Proof. Part (i). Continuous spectrum. Fix some A = -fj, 2m < and consider 

(3.13) B^ = AU = -fi 2m U in R N . 

Using Fourier Transform, it is not difficult to see that (13 . 1 3[) has a radial solution with 
the following asymptotic behaviour: as \y\ — > oo, 

(3.14) U A (y) ~ cos(ji\y\ + a ) {a eR). 

For m = 1, U\(y) is given by Bessel's function J u (/i\y\), with v = ^j^, [28, § 23]. Here, 
U A £ H 2m (R N ) (and £ H m (R N ), £ L 2 (R N )). Hence, (B 1 - A/)" 1 is not bounded for any 
A < 0, from whence G o~ c (Bi) by closure. 
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Discrete spectrum. We first discuss a general representation of all the eigenfunc- 
tions, which we will use later on. To this end, we introduce the polar coordinates x = (r, a) 
in R N , so that 

(3.15) A = A r + ^A CT , where A r = D 2 r + Z=± D r 

and A CT is the Laplace-Beltrami operator on the unit sphere S^ -1 in M. N , which is a 
regular operator with discrete spectrum in L 2 (S' Ar ~ 1 ) (each eigenvalue repeated as many 
times as its multiplicity) 

(3.16) cr(A CT ) = {-fi k = -k{k + N - 2) < 0, k > 0}. 

A CT has an orthonormal, complete, and closed subset {/fc(cr)} of eigenf unctions being k-th 
order homogeneous harmonic polynomials restricted to S' 7V_1 . 
Consider the eigenvalue problem for the operator (I3.2p . 

Bl U= (-l) m+1 (A r + £ A a ) m U + U 

(3 ' 17) = (-i) m+1 EZo C LK{^ ^T' l u + j^u = au, 

where C l m are binomial coefficients. Performing in (I3.17P separation of variables 

(3.18) U(r,a) = <j>(r)f k (<r) 
yields the purely radial eigenvalue problem for 4>(r), 

(3.19) (-l)^ 1 YZ, C l m (-Vk) m - l K <f) + = A0. 

Notice that the first operator in (13.191) contains singular potentials up to the leading 
singularity in the last term 

• •• — ^ 4>, with a constant c<i m = > 0. 

This and other singular terms have the right sign to guarantee that this operator in 
negative and coercive, and, of course, this essentially affects the final result. 

The positive lineal C + = Span{[/ , U x , ...} is infinite-dimensional. The rest of 
the results are standard in elliptic theory; see [2]. In particular, the existence of a countable 
set of eigenvalues follows from the fact that the subspace on which the quadratic form for 
the operator (13.61) . 

(3.20) Q 1 {v) = (v,B 1 v), 

is positive, is infinite-dimensional. Indeed, let £ C^°(R) be the cut-off function, 

1, < t < 1, 
0, t > 2, 



x(t) 



and set 

(3.21) Xa ,b(t) = x(t-b)[l -x(t-a)], 0<a<b =^ supp x ,6 e [a, b + 2]. 

We next chose a test sequence for Rayleigh quotion [2] in the following way: 

u a ,b(\x\) = |2r| -— ^Xo.bGzl). 
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Then we have 

Ql(u a ,b) = (u a ,b, BiM ,&) 
= (Ua,b, ( i+jfp - ]fpr)Wa,&) + ("0,6, [~ (~ A) m , Xa,&] 1 37 1 ) , 

where [-, •] is the commutator. Here, we have used the fact that (q.v. Section HJ) 



JV-2m 



-(-A) m \x\-— = j^\x\ 

so that the function [— (— A) m , x a ,b] m (13.22j) has support restricted to outer and inner 
layers according to (13.211) . Obviously, the first term on the right-hand side in (13.221) tends 
to +00 as b — > +00 since c > cn(m) and 

fi 1^ , 1 M _2m |x| _2JV 2 2m dx ~ f°° r N ~ 1 r - N dr = +00. 

In contrast, the second term on the right-hand side in (I3.22p remains bounded, since 

supp [-(-A) m , X a,b] C {a < \x\ < a + 1} U {b + 1 < |x| < b + 2}. 

As the result, Qi(w ,&) +°° as b — > +00, and there is a 6(a) such that Qi(w ffl ,6(o)) > 
and supp w aj fc G {a < |x| < 6+2}. Now choosing a sequence {a«} such that > 6(ai_i)+2 
for z = 1,2,..., we obtain infinite sequence {w 0< ,&(ai)} suc h that Qi(« ai ,6(ai)) > and 
supp u ai! b( ai ) H supp u a .b{ a ) = for « 7^ j. So there is an infinite dimensional subspace 
on which the quadratic form for the operator (13. 6p is positive, and consequently we have 
infinitely many positive eigenvalues of Bj. 

Part (ii). Exponential decay follows from separation of variables and ODE techniques. 
Part (iii). The positivity (13.121) for m = 1 is Jentzsch's classic theorem (1912). □ 

3.4. Nonexistence theorem. Thus, we consider the Cauchy-Dirichlet problem for ( 11. 12ft 
with general data ( 12.41) . Without loss of generality, we state the following nonexistence 
(divergence) result, where, for convenience, we slightly change the argument of the proof. 
This underlines extra features of the above positive lineal of the limit operator (13. 6p . 

(3.23) £+ = Span{[/ ,t/i,...}, 

as a set of all finite linear combinations of the given eigenfunctions. 
Theorem 3.1. Let 

(3.24) u 0£ {ey) -> v {y) as £ — >■ m Lf oc (R JV ), 
and let, in the metric of L 2 (R N ), 

(3.25) v is not orthogonal to C + . 
Then (TO} holds. 

The proof repeats the arguments of the proof of Theorem 12.11 In particular, assuming 
the continuity and strict positivity at the origin of the data that are independent of e, 

(3.26) uo e (x) = uo(x) in Q, and Uq G C(Bi), wo(0) = So > 0, 
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we obtain from (I3.24[) 



(3.27) v (y) = S = u (0) > 0. 

Then the condition (I3.25P (more precisely, see (i) below) is obviously valid for m — 1 in 
view of the positivity (13.121) . For m > 1, checking (13. 25j) is not straightforward since all 
the eigenf unctions are of changing sign. Anyway, we expect that (13.251) is always valid for 
constant data (13.271) . It is curious that, for m > 2, this is an open problem. 

3.5. On special solutions with fixed nodal sets. Using the separation ( 13.171) . we look 
for special solutions 

(3.28) u(x,t) = v(r,t)f k (a), r = \x\, 

of the original singular equation (II. 6p with Q — B^. These solutions have fixed (time- 
independent) nodal sets that actually changes the Hardy constant since now effectively 
belongs to an artificial boundary, on which u = 0. Indeed, formally substituting ( 13.281) 
into (jl.6p yields the following radial equation for v(r,t): 

(3.29) v t = (-l^+i^cU-^r-'Ki^Lwv) + ^v. 
It follows that the last most singular at y = term is now 

(3.30) + 

so that the problem becomes subcritical provided that 

(3.31) c - (/i fc ) m < Oh, i.e., for c < c H + (^) m 

Thus, for any c > ch, there exist global solutions ( 13.281) of ( 11.61) with k 1. 

Note that this by no means undermines the phenomenon of the strong instability in 
the singular equation (II. 6p in the supercritical range (ll.9p . Indeed, it is not clear which 
solutions (I3.28p . (I3.3ip do withstand the e-regularization without blow-up as e —> 0. We 
expect that, without the assumption Q = B\ or without special symmetries of Q sup- 
porting the fixed nodal set of special solutions under consideration, no solutions actually 
do. However, in the maximal generality, the proof is not a part of our business here. In 
Section [5j we present the proof of blow-up for some stationary data that are supposed to 
be the most resistive to regular e-approximations, but actually are not. 

4. Constructing data with doubly oscillatory ±oo limit 

In this section, for the reason of performing some rather involved calculus, we take the 
unit ball, Q = Bi, and restrict to radial solutions only, where r = \x\ & [0, 1). 
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4.1. Singular stationary solution (SSS). Consider the radial stationary equation 

(4.1) Bo U = -(-A r ) m U+^U = for rG (0,1). 
As usual for Euler's-type ODEs, looking for solutions of (14.11) 

(4.2) U{r) = r 7 , 

we obtain the following characteristic equation for 7 G C: 

G{i) = G>(7) + c = 0, where 

' U ' G.( 7 ) = R [7 - 2(fc - 1)](7 + iV - 2*). 

fc=i 

Then (jl.lOp means that (see e.g., [8]) 

(4.4) c = CH = _ a (-^)> , 

so the best Hardy constant ch is such that the function (14. 2 j) with the exponent 
(4-5) 7 = 7m = 

is the exact weak solution of the homogeneous equation Bo^ = in R . Thus, for c = ch, 
the characteristic equation has the double root (14.51) . It generates two L 2 -solutions, 

— N-2m — N — 2m 

(4.6) U m (r) = r 2 lnr and U m+ \{r)—r 



2 



which are ordered relative to the growth rate as r — > 0. Other 2m — 2 characteristic roots 
of ( 14. 3 ft are real or complex. Complex roots can occur for m > 3. One can see from the 
structure of the characteristic polynomial (14. 3 p that, for iV ^> 2m and m even, there exist 
precisely two more real roots 7 m _i > 2(m — 1) and 7 m _i < 2 — N. On the other hand, 
for N ^> 2m and m odd no more real roots exist; see [8]. 

It follows from (14. 3p that in the supercritical range c > ch, the characteristic polynomial 
admits two complex roots, 

(4.7) 7m± = -^^±irf, where d = 0(^=c^) > 0. 
The corresponding solutions are oscillatory near the origin, e.g., 

N — 2m 

(4.8) U m (r) = r 2 C os(<ilnr + a ) (a G R). 

By a weafc 5'S'5' of ( 14.11) denoted by U*(r), we mean a solution that exhibits the oscillatory 
behaviour such as (14.81) as r — > 0. Note that 



(4.9) [/* G L p (fi) for any p < 



2N 



N-2m ' 
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4.2. e-approximation diverges to ±00. In the radial parabolic approximating problem 
for equation (II. 121) . we first take initial data Uo(r) that are not continuous at r = and 
exhibit an oscillatory behaviour. It is convenient to assume the behaviour as in (14.81) . i.e., 

(4.10) Uo(r) = r 2 cos(<ilnr) for r > small (or, simply, Uo(r) = U*(r)), 

though other types of oscillations will also do. 

Let us explain the origin of such huge oscillatory properties of the approximation 
u £ (x,t). Using the same idea, we apply the estimates from Section [2] that are sufficient. 
The simpler case m = 1, where U (y) > (q.v. (13.121) ) was studied in [9J, so here m > 2. 

Thus, we begin with formal estimates of the first Fourier coefficient for the maximal 
positive eigenvalue as in (12 . 1 5[) (here r = \x\ > 0), 

(4.11) X e « ^ and ipo(x) ~ £~^^7 (f ) ~ e"f e ~ amA o /2m I cos (b m ^) for r > e, 

where /i m = | (— a m + ib m ) is the root of (— l) m+1 /i 2m = with the maximal Re/x m = 
— a m < 0. This yields, in the first rough, but sufficient, approximation, the coefficient 

1 1 

(4.12) cl= (u ,r ) ^e-T j r^uoi^Uo^) dr = ef f z N - 1 u {ez)U {z)dz. 



Substituting next data (I4.10p yields the following estimate of the first Fourier coefficient: 

+00 

(4.13) c £ ~ e m J z E±2 ^U (z) cos(d\nz- D) dz = e m I d (D), 



where D = — (ilne-^+ooase-^Oisa shifting argument in cos(-), which changes the 
sign of the function in the integral. Recall that Uq(z) is exponentially small by (14.111) as 
z — > 00, so the integral in (I4.13P fast converges. Here, we assume that both integrals 

+OO Ar+2 2 N+2 2 

f z 2 U (z) cos(d\nz) dz and f z 2 U (z) sin(<iln z) dz 


do not vanish simultaneously (otherwise, we change d and/or the exponent — N ~^ m or any 
other parameters to get necessary non-zero values). Then, in view of analytic dependence 
on parameters in (14.131) . for any fixed d > 0, there exists sequences {e^ } — > such that 

c k > (£fc) m £>o and c fc < — (£/T) m <$o for any k = 0,1,2,... with some constant 5$ > 0. 
Therefore, by fl4TTll . for e = — > 0, 

(4.14) u £ (x,t ) ~ c%e m -^ I d (d\ne)e~ amA ° 2m i cos(6 m f) e A ° & (r = \x\), 

where we may assume that cos(-) does not violate the above sigh restrictions on c £ Q . 
Similarly, in view of the oscillatory functions cos (b m r ) in (I4.14p . there exist sequences 
{i^} — > such that 4>q(x) stands positive and negative respectively. 

Thus, we observe a doubly oscillatory singular limit to ±00 as e — > for the regularizing 
solution ( 14.14p . Imposing necessary elementary assumptions on these four sequences {e^ } 
and {e^ }, one always can choose some "intermediate" sequences {e^ } — > along which 
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the limits ±00 are guaranteed. Thus, by (I4.14p . we observe that, for any r = \x\ > 0, in 
the sense of the first Fourier coefficient, 

(4.15) limsup £ ^ u e (x, t) = +00 and liminf e ^ u E (x, t) = — 00 
In particular, we also easily fix the following weaker ±oo-divergence: 

(4.16) limsup e _> sup,,. u e (x,t) = +00 and liminf e ^ infz w £ (x,£) = —00 

It is interesting that, by oscillatory behaviour of the Fourier coefficient Cq in (14.131) . the 
non-uniform divergence (14.151) holds even at the origin x = 0. 

On the other hand, if for given data Uq, the coefficient Cq is not oscillatory at all, e.g., 

Cq > c £ > for all small e > 0, 

the oscillatory divergence (14.161) or (I4.15P for any x 7^ holds due to the sign-changing 
behaviour of U (y) for y ^> 1 given in (14.111) . 

5. An analytic stationary profile blows up as e — > 

Let us present a different simple example showing another unstability feature of the 
^-approximation of (11.61) . which we now consider in R N x (0,T), i.e., for convenience, we 
pose the Cauchy problem. We choose c > ch such that the characteristic equation (14.31) 
has the root 7 = 2m, so that (14.11) has the analytic stationary solution 

(5.1) U (x) = \x\ 2m in R N . 

Consider the regularized equation (11.121) with the stationary data (15. ip . Choosing, for 
convenience, the new variable 

v = u t 

and differentiating (11.121) in t, we obtain for v the same equation (I1.12p with data 
v 0£ (x) = B £ \x\ 2m = -(-A) m \x\ 2m + g2m+ ^ |2m \x\ 2m 

( 5 - 2 ) — ( 1 1 \ _ ce 2m 

~ \s 2m + \x\ 2m \x\ 2m ) e 2m +\x\ 2m ' 

After scaling (13.11) . we arrive at the equation (13.51) with e-independent initial data 
(5.3) v (y) = -j^ < in R N . 

Therefore, (I3.25P guarantees blow-up as e — > 0. As usual, m = 1 is easy by (13.121) . 

Thus, in general, even classical stationary solutions of the singular equation (I1.6P do 
not stand e-regularization, to say nothing about other more arbitrary (non-steady) data. 

6. Parabolic extensions: time-dependent and nonlinear potentials 

6.1. Time-dependent potentials: towards non self-adjoint spectral theory and 
a new Hardy constant. New spectral phenomena appear when the potential in (12.11) 
depends on the time- variable in an e-scaling manner, e.g., as in equation with the following 
singular critical potential (cf. (II. 6p ): 

(6.1) u t = -(-A) m u + JT ^u. 
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Then the natural regularized equation reads 

(6.2) u £ : u t = —(-A) m u + £2m+ ° Mx?m u. 

As usual, we treat a general regularized equation with such potentials, 

(6.3) u £ : Ut = B £ u = -(-A) m u+^q(^ 7 L)u in R N xR + . 

For convenience and as a different example, we now ignore boundary conditions and 
consider the Cauchy problem with data u 0£ G L 2 {R n ) H L°°(M jv ). Scaling (13.11) yields for 

(6.4) u £ (x,t)=v £ (z,^) 
the following rescaled equation: 

(6.5) v £ : v T = ~(-A) m v + q(y,r)v, v(y,0) = v 0£ (y)=u 0£ (ey). 

Since by ( 16. 4ft the behaviour of {u £ } as e — > is equivalent to the behaviour of {v £ (y, r)} 
as r — > +oo, we perform extra scaling and introduce the new standard similarity variable 
from the heat kernel of the operator D t + (— A) m : 

(6.6) z = (1+r )i/ 2m , s = ln(l + t) -> +oo. 
Then u> e (z, s) = v £ (y,r) solves the equation 

(6.7) w £ : w s = — (— A) m w + -2 ■ Vw + e s q(ze^, e s — l) w. 

In order to approach a fixed regularized potential (as in (I3.5P ). we assume that 

(6.8) e s q( ze?™,e s — l) — ► 1+ |^2m as s - +°° uniformly 
and sufficiently fast. For instance, (16.81) holds for the potential in (16. 2p since 

(6-9) ^2m ?( j) ^2m) £ 2m | x 1 2m '* lilj 1 ^) l + T+|?/| 2m ' 

Thus, passing to the limit in ( 16.71) with the assumption ( 16. 8p . we obtain the equation 
(6.10) w s = -{-A) m w + ^ z ■ Vw + j^w. 

Here the limit operator 

(6-H) B 1 = -(-Ar + ±z-V + J ^I 

admits Friedrichs self-adjoint extension for the second-order case m = 1 only, where it 
possesses a symmetric representation of the form 

(6.12) B 1 = iV-(pV) + TT ^I in L 2 p (R N ), with p(z)=e i t. 

For m > 1, (16.111) is not self-adjoint and is not symmetric in any weighted L 2 -spaces. 
Consider the linear operator with the same principal differential part, 

(6-13) B = -(-A)™ + £z.V+&I. 

It is known [6j that, for any m > 1, operator ( 16.131) is naturally defined in the weighted 
space LJ(R*), with 

(6.14) p( y )=e a ^ a , with a = ^rr, 
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where a > is a sufficiently small constant. The domain of B is the corresponding Sobolev 
space H% m (R N ), so that B : H* m (R N ) -> L 2 p m (R N ) is bounded with discrete spectrum 

(6.15) (T (B) = {-^ 1 = 0,1,2,...}. 

In particular, both B and B* then have discrete spectra, compact resolvents, and complete 
sets of eigenfunctions. Some of these results can be extended to the operator (16. lip , for 
which spectral theory deserves further development. 

We now concentrate on particular spectral properties related to the nonexistence. 
Consider first the self-adjoint case m = 1 . Then, by classic theory [21 [20], there 
exists a monotone branch Ao = Ao(c) of the first simple real eigenvalue of Bi (here m — 1): 

(6.16) A (0) = — and A (c) — > +oo as c — > +oo. 

We then define the new "Hardy constant" for the c-family of operators (16.111) as follows: 

(6.17) c H : A (c) = 0ea(B 1 ). 

Then, for any c > ch, there exists A (c) > that leads to generic blow-up of approxi- 
mations {u £ } provided that v is not orthogonal to U . Analogously, further increasing 
c > c H leads to appearance of more branches of positive eigenvalues Ax(c), A K (c), so 
that we can have the positive lineal C + = Span {U , Uk} of arbitrarily large, but finite 
dimension K = K(c) for c ^> 1. 

Thus, taking into account the first unstable mode only, we have from scalings (16.41) 
and (16.61) that, under the "non-orthogonality assumption" to Uq G £+, the following 
divergence rate is achieved: for an arbitrarily small fixed t > (here again m — 1), 

(6.18) u E (x,t) ~ (^) A(,(c) [/ (^) as e - 0. 

Similar to Theorem 12.11 the behaviour (16.181) is easy to characterize in terms of L P (M. N )- 
divergence for any p > 1, as well as in other metrics. 

On the contrary, for c < Ch, the family {u e = v e = w e } is uniformly bounded meaning 
existence of a bounded solution (of (16. ip . say). 

Consider m > 2, where we can guarantee for sure less. Since in (16. lip the last term 
1+ l^2m I serves as a compact perturbation of B in the equivalent integral representation, 
we still have a discrete spectrum and a countable family of continuous eigenvalue branches 
for all c > 0. For c = 0, there holds: 

B l = B ~ 2^ J ' 

so by perturbation methods [20], we find the local branch as in ( 16. 161) of the first real 
simple eigenvalue Ao(c) of Bi for all sufficiently small c > 0. Further extension of the 
branch in this non self-adjoint case is not guaranteed to be strictly monotone, and the 
eigenvalue does not necessarily remain always real. Then we are assumed to deal with 
the first eigenvalue having the maximal real part (q.v. (16.171) ). 



(6.19) 



Re A (c) > for c> c H . 
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This definition is based on a natural property that ReA (c) growths unboundently as 
c — ► +00, though passage to the limit c — > +00 in the family (16. lip is quite tricky and 
not studied here. Then under the assumption ( 16.191) (or ( 16. 17ft in the real case), we still 
obtain the divergence result of the type (16.18P for any non-orthogonal initial data now in 
the sense that {vq, Uq) 7^ 0, where Uq is the corresponding adjoint eigenfunction of B^. 

Thus, for all m > 1, according to our previous results, the divergence of the family 
{u £ } as e — > will be guaranteed under the hypothesis of the type (13.251) . where C + is 
the positive lineal of the limit operator (16.111) . Indeed, deeper spectral analysis of the non 
self-adjoint operators (16. lip are quite in demand here. 

6.2. On nonlinear parabolic equations. As a next natural extension, we consider the 
following semilinear parabolic equation with the same singular potential as in (11. 6p : 

(6.20) u t = -{-A) m u + W here p > 1. 

Here, the nonlinearity is taken in the form \u\ p for convenience of further simple calculus 
in the case m > 2, when the Maximum Principle fails. 

The condition p > 1 involves the extra high instability of the evolution driving by the 
reaction-diffusion PDE (I6.20p . since now it described blow-up of solutions in finite time. 
For m = 1 and u > 0, the nonexistence results for (16.201) were obtained in Brezis-Cabre 
[3J; on full history and extensions see [HI [22], [3 P- 333], and [23 p. 108,267]. For 
solutions of changing sign, the nonexistence results for ( I6.20p were absent even for m = 1. 

Let us see how the nonlinearity \u\ p affects the regularized solutions satisfying, as usual, 

(6.21) u £ : u t = -{-A) m u+ £2m ^ x?m \u\v. 
Then scaling ( 13.11) yields 

(6.22) v £ : v t = -(-A) m v+ TT ^\v\P in R N x R + , 

and for simplicity we assume that data are constant as in (13.271) . 

Case m = 1. Then the proof of blow-up of (16.221) is straightforward, since by (13.271) . 

(6.23) v t > in l ff x K + 

via the Maximum Principle. Therefore, the solution v(y, r) is strictly monotone increasing 
in t and is not bounded above (since otherwise, as for gradient systems, v(r) would 
stabilize to a bounded steady solution, which is nonexistent); see details in [7,, p. 319]. 
Denoting by T* the blow-up time for ( 16. 22ft . we take into account the fact that due to 
(16.231) blow-up is complete, i.e., in the natural sense of the proper minimal extension, 

v(y, t) = +00 for r > T*. 

We then conclude that, for all small e > 0, 

(6.24) u £ (x,t) blows up completely in time T £ < 2e 2m T*. 

This ends the proof of divergence as e — > of {u £ } under the assumption (I3.26P for m = 1. 
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Case m > 2. A similar scheme applies with some changes since f!6.23|) is no longer true. 
Finite time blow-up in (16.221) . (13.271) is then proved by the eigenfunction method that 
shows blow-up of the first Fourier coefficient 

where tpi(z) is a specially constructed nonnegative cut-off function satisfying the following 
elliptic inequality: there exists a constant Ai > such that 

(6.25) |A">i| < AiV>i in R N . 

See a convenient adaptation of Kaplan's eigenfunction method (1963) to higher-order par- 
abolic operators in [TUl § 5]. It is crucial that the resulting ordinary differential inequality 
for E(r) implies not only the fact of finite-time blow-up at some T* but also that it be- 
comes complete (possibly, at some moment T c > X 1 *). An alternative proof of blow-up in 
(16.221) by the nonlinear capacity method can be performed as in [251 § 29]. This leads to 
the nonexistence conclusion (16.241) . Here, similar to (I2.15p . we control complete blow-up 
of the "first Fourier coefficient" relative to a solution ip% of the elliptic inequality (16. 25ft . 
The above scheme on nonexistence applies to the truly quasilinear parabolic model 

(6.26) u £ : ut = -(-A) m (\u\ n u) + g2m+ 1 | , |2m \u\ n+1 , where n > 0. 

7. On nonexistence for other linear and semilinear PDEs 
7.1. Schrodinger-type equations. A typical linear example is 

(7.1) u £ : iut = B £ u = -(-A) m u + g2m+ % |2m u. 
The principle difference is that all u £ satisfy the conservation law 

(7.2) K(*)||l» = \\uoe\\v> for all t > 0. 

Nevertheless, the nonexistence conclusion for c > ch is derived in a similar manner, 
though, in view of (17. 2p . the key feature of divergence is that the family {u £ (x,t)} gets 
extremely oscillatory as e — > in the variable t > 0. Indeed, now, loosely speaking, for 
positive eigenvalues Aj of B 1; 

(7.3) X £ j ~ — i for small e > 0. 

Finally, instead of (12.151) . in view of the oscillatory behaviour with purely imaginary 
eigenvalues (|7.3|) . assuming the strict inequality A > A 1 (cf. ( I3.10p ). we suggest the 
following characterization of the actual divergence of solution (of course this is just a 
formal illustration that is not that informative, which any isolated mode admits): 

(7.4) e iA §*M £ (x,t) - cg^o(f) ^ as e -> weakly in Ll c (Q x (0,5)). 

Thus, in terms of u £ , we get rather weak blow-up. However, the time-derivative (u £ ) t is 
then truly blows up as e — > according to 

K)t(^*)~iLe- iA S*cg^(f). 
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The real parts of the solutions, u £ = U £ + iV £ , satisfy a 4mth-order hyperbolic equation, 

U e : U tt = -B 2 e U, 

to which a general scheme of the analysis, with some changes, can be also applied. 

As usual, a similar "nonexistence" analysis of the corresponding nonlinear Schrddinger 
equation (NLS) 



(7.5) u £ : iut = B £ u = -(-A) m u + £2m+ 1 |a;|2m \u\ p x u (p > 1) 
demands extra results on blow-up in the corresponding rescaled equation 

(7.6) v £ : i v T = B lV = -{-A) m v + TT ^ r \v\*~ x v. 

For m = 1, blow-up for such NLSEs in some parameter ranges of p and iV is a classic 



matter of the theory and application to nonlinear optics, plasma physics, and others; see 
references to key papers and monographs on this subject in Merle-Raphael [23] and in 
survey [12]. The higher-order case m > 2 is less developed though some techniques can 
be applied. Notice that we are interested in quite particular setting with constant data 
as in (13.271) . for which the proof of blow-up seems not very difficult. 

Thus, we will get typical "divergence" of the family {u e } as e — > 0, provided that the 
suitable rescaled solution of (17. 6p blows up in finite time. Note that finite L 2 -energy blow- 
up solutions of NLSEs admit a natural extension beyond blow-up time, i.e., for t > T 
in view of the conservation (17.21) (possibly, this is not the case for data (13.271) of infinite 
energy). Then the divergence of {u £ } as e — > might be not that lethal as in other 
examples, so again a further study of such singularity phenomena is necessary. 

7.2. Hyperbolic equations. A key regularized linear model is constructed analogously, 

(7.7) u £ : u tt = B £ u = -(-A) m u + e2m +\ x \2 m u, 

with initial data u , u\. Taking t = e m r in (13. ip yields the rescaled equation 

(7.8) v £ : v TT = B 1 v = -(-A) m v + TT ^v. 

The separation of variables generates the eigenvalue problem (12. 9p . where A i— > A 2 . Hence, 
the positive lineal of Bi continues to play the key role for the nonexistence as e — > 0, so 
under similar non-orthogonality assumptions, we arrive at the divergence result such as 
( 12.151) . with v/Af in the first multiplier. 

The corresponding nonlinear hyperbolic equation, 

(7.9) u £ : u tt = B £ u = -(-A) m u + £2m+ 1 |i p m \u\ p , 
leads to the rescaled PDE 

(7.10) v £ : v TT = B 1 v = -(-A) m v + JT ^\v\P. 

For data (I3.27p . the necessary blow-up results for (I7.10p can be obtained as in [251 Ch. 5-7], 
so that the nonexistence-divergence conclusions as e — > persist in a similar way. 
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